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Preface 


The  purpose  of  this  project  was  to  determine  the 
position  errors  of  the  angle-of-attack  (AOA)  sensors  on 
aircraft  using  state  estimation  with  flight  test  data. 

Aircraft  from  the  USAF  Test  Pilot  School  (TPS)  were  used  to 
obtain  flight  test,  data,  and  Kalman  filtering  was  used  to 
process  the  data.  The  results  of  this  project  are 
significant  to  future  flight  test  projects  where  an  accurate 
AOA  measurement  is  required. 

Aircraft  AOA  position  errors  are  caused  by  aerodynamic 
factors  such  as  local  flow  and  upwash.  The  first  step  in 
finding  those  errors  was  to  determine  the  equations  for 
calculating  the  true  AOA  from  other  available  flight  test 
parameters.  Since  the  inputs  to  those  equations  were  from 
instrumentation  on  flight  test  aircraft,  they  were  noise 
corrupted  and  had  to  be  filtered.  I  used  state  estimation 
in  a  Kalman  filter  program  to  calculate  an  "optimal"  true  AOA. 
The  data  were  obtained  from  flights  in  a  T-38A  Talon,  a 
two-seat,  supersonic  trainer  modified  with  an  instrumented 
Vought  yaw  and  pitch  system  noseboom.  The  position  errors 
calculated  in  this  report  are  only  good  for  that,  aircraft  and 
nose  boom  configuration.  However,  the  methods  used  are 
applicable  to  all  properly  instrumented  aircraft. 

I  would  like  to  thank  my  thesis  advisors.  Major  (Dr.) 
James  T.  Silverthorn  of  the  USAF  TPS  and  Dr.  Robert  A.  Calico 
of  AFIT,  for  their  help  in  this  project.  I  would  also  like  to 


thank  the  test  pilots  I  flew  with  on  the  data  flights.  Major 
Philip  B.  Arnold  and  Captain  David  J.  Eichhorn  of  the  USAP 
TPS.  But  most  of  all  I  would  like  to  thank  my  wife,  Diana, 


for  her  help  and  understanding  over  the  last  two  years  of  AFIT 
and  TPS. 
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Abstract 


This  project  determined  the  position  errors  of  an 
aircraft's  angle-of-attack  (AOA)  sensor  using  state 
estimation  with  flight  test  data.  The  position  errors  were 
caused  by  local  flow  and  upwash  and  were  found  to  be  a 
function  of  AOA  and  Mach  number.  The  test  aircraft  used  in 
this  project  was  a  T-38A  Talon  supersonic  trainer  from  the 
USAF  Test  Pilot  School  configured  with  a  Vought  yaw  and  pitch 
system  noseboom  and  an  internal  Aydin-Vector  data  acquisition 
system  (DAS). 

The  position  errors  were  found  by  calculating  the  true 
AOA  using  equations  of  motion  and  DAS  parameters.  The  data 
from  the  DAS  were  noise  corrupted  and  had  to  be  filtered. 

This  was  accomplished  using  state  estimation  in  a  Kalman 
filter.  The  estimated  AOA  was  compared  to  the  measured  AOA 
from  the  noseboom  sensor  to  obtain  the  position  error. 

Accurate  position  errors  were  obtained,  even  in  dynamic 
maneuvers.  The  method  was  accurate  enough  to  identify  a  hyster¬ 
esis  error  in  the  T-38A's  AOA  sensor  of  +/-  0.5  degrees,  which 
was  confirmed  by  ground  calibration.  This  method  should  be 
considered  in  future  AOA  error  testing. 
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USE  OF  STATE  ESTIMATION  TO  CALCULATE 
ANGLE-OF -ATTACK  POSITION  ERROR 
FROM  FLIGHT  TEST  DATA 


I.  Introduction 


Problem 

Angle-of-attack  (AOA)  is  a  primary  parameter  of 
performance  and  stabi li ty-and-control  in  flight  test. 
Unfortunately,  the  AOA  measured  by  the  aircraft  sensors  has  a 
position  error  caused  by  the  aerodynamic  influence  of  the 
aircraft  body.  The  first  source  of  this  position  error  is 
local  flow  about  the  AOA  sensor  caused  by  aerodynamic 
interference,  boundary  layer  effects,  and  shock  interaction. 
The  second  source  of  the  position  error  is  upwash  from 
aircraft  components  such  as  the  fuselage  and  wing.  The 
accurate  determination  of  AOA  position  error  is  a  significant 
problem  in  flight  test  (1:7). 


Background 

The  magnitude  of  the  AOA  position  error  is  evident  with 
the  USAF/CAL  variable  stabili  ty  NT-33  airplane,  a  jet  trainer  v-| 

used  in  the  USAF  Test  Pilot  School  (TPS)  curriculum.  The 
NT-33  has  a  fuselage-mounted  AOA  vane  which  is  subject  to 

large  flow  and  upwash  effects.  Figure  1  shows  the  NT-33  AOA  . 

(and  sideslip,  which  has  similar  errors)  position  error  /- 
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Wind  tunnel  calibration  is  a  commonly  used  method  of 
determining  AOA  position  error  at  the  Air  Force  Flight  Test 
Center  (AFFTC).  An  entire  noseboom  instrumentation  unit 
(NBIU)  can  be  installed  in  a  wind  tunnel  and  tested  over  a 
range  of  conditions.  One  of  the  AFFTC  standard  NBIUs,  a 
Conrac  adapter  with  a  Rosemount  Model  852G  pitot-static  probe, 
was  tested  in  the  NASA/Ames  Research  Center  wind  tunnels  in 
1973  (3).  The  AOA  position  error  was  found  to  be  small,  less 
than  7.5%  (3:38).  The  wind  tunnel  test  showed  AOA  position 
error  to  be  a  function  of  Mach  number  and  sideslip  angle. 
Reynold's  number  effects  were  not  discovered. 

Very  little  flight  testing  has  been  accomplished  to 
determine  AOA  position  error.  One  technique  that  has  been 
used  is  to  mount  flight  path  accelerometers  on  the  test 
aircraft  and  fly  1  g,  wings  level  stable  points  over  a  range 
of  Mach  numbers,  sideslips,  and  aircraft  weights.  True  AOA 
(«)  is  determined  from  the  equation: 

a  =  0  -  y  ( 1 ) 

where  e  is  pitch  angle  in  wings  level  flight  and  y  is  the 
flight  path  angle.  The  T-46  jet  trainer  Combined  Test  Force 
is  planning  to  use  this  technique  to  calibrate  their  AOA 
sensors  when  flight  testing  begins  in  October  1985.  They  plan 
to  eliminate  noise  in  the  data  by  using  a  2  hz  low  bypass 
Butterworth  filter.  Unfortunately,  not  all  aircraft  can  be 
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equipped  with  flight  path  accelerometers  due  to  the  size  and 
required  cost-  Furthermore,  the  range  of  AOA  that  is 
attainable  at  a  particular  Mach  number  is  very  limited  for 
straight  and  level  flight,  since  altitude  is  the  only  variable 
that  can  be  adjusted.  As  an  example,  the  T-38A,  at  Mach  0.83, 
flies  at  +  2.5  degrees  AOA  at  25,000  feet  and  at  +  1  degrees 
AOA  at  15,000  feet.  A  technique  to  obtain  AOA  position  error 
during  dynamic  maneuvers  is  required. 

A  new  flight  test  method  of  determining  AOA  position 
error  is  through  the  use  of  MMLE3,  a  modified  maximum 
likelihood  estimation  program  (4)  MMLE3  uses  the  aircraft 
mathematical  model  with  estimated  stability  and  control  (S&C) 
derivatives.  Flight  test  maneuvers  such  as  elevator  doublets 
are  flown,  and  MMLE3  tries  to  match  the  time  history  of  the 
maneuver  with  the  time  history  of  the  math  model  by  changing 
the  estimated  S&C  derivatives.  MMLE3  also  calculates  an  AOA 
position  error  factor  for  the  maneuver  (4:3).  MMLE3  is  not 
extremely  accurate  and  requires  numerous  flight  test  maneuvers 
to  increase  its  accuracy.  An  easier  and  more  accurate 
technique  is  needed  to  calibrate  AOA  sensors. 

Scope 

The  purpose  of  this  project  was  to  determine  the 
position  errors  on  the  AOA  sensors  of  aircraft  using  flight 
test  data  available  from  standard  data  acquisition  systems 
(DAS).  Initially,  AOA  position  error  was  determined  using 
deterministic  equations  from  straight  and  level  flight. 
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Problems  with  this  technique  suggested  a  more  general 
approach.  State  estimation  in  the  form  of  Kalman  filtering 
was  used  to  filter  out  noise  on  the  flight  test  data  and 
calculate  an  "optimal"  true  AOA  during  dynamic  maneuvers. 
This  true  AOA  was  compared  to  the  measured  AOA  to  determine 
the  position  errors.  USAF  TPS  T-38A  aircraft  were  used  to 
collect  data  and  the  AOA  position  errors  are  valid  for  those 
aircraft.  However,  the  technique  will  work  for  any  properly 
instrumented  aircraft. 

Objectives 

The  objectives  of  this  project  were  to: 

(1)  Determine  the  equations  necessary  to  calculate  the  true 
AOA  from  flight  test  data. 

(2)  Use  state  estimation  (Kalman  filtering)  to  filter  noise 
from  the  flight  test  data  and  calculate  an  "optimal"  tr 
AOA. 

(3)  Collect  the  flight  test  data  needed  to  compute  the 
"optimal"  true  AOA. 

(4)  Calculate  the  AOA  position  error  correction  factors  for 


the  test  aircraft. 
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II.  Angle-of-Attack  Equations 


AOA  Correction  Factor 


The  AOA  position  error  correction  factor,  K  ,  is 

a 

calculated  from  the  equation  (4:3): 


q  +  ‘"h 


P  (ya) 


where  a  is  the  true  AOA  of  the  aircraft  and  a  is  the  AOA 
c  m 

measured  by  the  aircraft's  sensor. 


The  term : 


q  (xa  +  Leg } 


corrects  the  measured  AOA  for  pitch  rate  (q)  effects.  The 
terms  x^  and  Acg  account  for  the  longitudinal  distance  from 
the  eg  to  the  AOA  sensor.  U  is  aircraft  true  airspeed. 

The  term: 


P  (Y«) 


corrects  the  AOA  for  roll  rate  (p)  effects.  The  term  is 
the  lateral  distance  from  the  aircraft  centerline  to  the  AOA 
sensor . 

All  flight  testing  for  this  project  was  accomplished 
wings  level.  Since  there  was  no  roll  rate,  equation  (4)  drops 
out  of  equation  (2).  Pitch  rate,  true  airspeed,  and  measured 
AOA  are  parameters  measured  by  the  aircraft  DAS.  Longitudinal 


distance  from  the  eg  to  the  AOA  sensor  is  a  function  of 
aircraft  fuel  weight  and  is  easily  calculated.  The  only 
remaining  unknown  is  true  AOA. 


Equations  of  Motion 

In  order  to  calculate  the  true  AOA,  an  equation  was 
needed  that  used  parameters  available  from  the  aircraft  DAS. 
Equation  (1)  showed  the  angular  relationship  between  AOA, 
pitch  angle,  and  flight  path  angle  in  wings  level  flight. 
Since  many  test  aircraft,  including  the  TPS  aircraft  used  in 
this  project,  do  not  have  flight  path  accelerometers,  flight 
path  angle  ( y )  must  be  calculated  by: 


sin  r 


(5) 


where  h  is  the  vertical  velocity  of  the  aircraft.  Vertical 
velocity  can  be  calculated  as  the  time  rate  of  change  of  the 
altitude  from  the  DAS. 

Another  equation  to  calculate  true  AOA  comes  from  the 
aircraft's  equations  of  motion  (5:3.21-3.51).  The  vector 
equation  for  applied  force  (F)  is: 


F 


m 


dV 

dt 


XYZ 


(6) 


which  applies  to  inertial  space.  Assuming  the  forces 
resulting  from  the  earth's  rotation  and  coriolis  effects  to  be 
negligible,  a  fixed  earth  axis  system  can  be  used  instead  of 
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inertial  space.  The  movement  of  a  vehicle  with  respect  to  a 
fixed  earth  axis  is  shown  in  Figure  2. 
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Figure  2.  Relationship  of  Fixed  Earth  Axis  (XYZ)  to  Vehicle 
Body  Axis  (xyz)  (5:3.22) 


The  vector  equation  for  the  time  rate  of  change  of 


velocity  from  one  axis  system  to  another  is: 


+  u  X  V, 


where  XYZ  is  the  fixed  earth  axis  and  xyz  is  the  aircraft  body 
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Equation  (6)  now  becomes: 


where  aircraft  velocity  (VT)  can  be  written  as: 

Vm  =  ul  +  Vj+Wk 
T  J 

and  aircraft  rotation  (u)  can  be  written: 

u  =  Pi  +  Qj  +  Rk 
Equation  (8)  now  becomes: 


Taking  the  cross  product  of  the  inner  term  and  expanding: 

F  =  m  CUT  +  V]-  +  WJ<  +  (QW  -  RV)T  -  (FW  -  RU)J  +  (PV  -  QU)k] 

Looking  at  only  the  z  component  of  force  gives: 

F  =  m  (W  +  PV  -  QU)  =  m  (a  ) 
z  z 

where  (W  +  PV  -  QU)  equals  the  normal  acceleration,  a  • 

z 

Assuming  that  the  aircraft  motion  consists  of  small 
deviations  from  an  initial  reference  condition,  the  above 


Q  =  Qo  +  q 

U  =  U  +  u 
O 

where  the  small  case  values  are  the  small  perturbations  from 
the  initial  values.  Assuming  the  aircraft  starts  from  wings 
level,  steady  straight  symmetrical  flight: 

W  =  P  =  V  =  Q  =  0 

o  o  o  o 

Equation  (11)  now  becomes: 


m  (w  +  pv  -  qU)  =  m  (a  ) 


(12) 


All  testing  during  this  project  was  done  wings  level,  so  roll 

rate  (p)  is  zero.  The  change  in  velocity  (u)  is  assumed  to  be 

small,  so  U  =  U  .  Dividing  each  side  by  m(U  )  gives: 

o  o 


(13) 


Assuming  small  AOA  gives  the  relationship: 


(14) 


which  can  be  substituted  into  equation  (13)  and  rearranged  to 


where  both  pitch  rate  and  normal  acceleration  are  measured  by 


the  aircraft  DAS.  Assuming  a  finite  time  interval  At,  a  =  Aa. 
An  iterative  equation  can  be  formed  where: 

+  ^  =  on  +  A  a  (  At )  (16) 

Computer  Program  AOAOPT 

A  FORTRAN  computer  program  was  designed  to  calculate 
true  AOA  using  both  the  angle  relationship  (equations  (1)  and 
(5))  and  the  iterative  relationship  (equations  (15)  and  (16)). 
This  program  is  called  AOAOPT  and  is  shown  in  Appendix  B.  The 
program  reads  the  required  flight  test  parameters  from  the 
DAS,  makes  necessary  pitot-static  corrections,  calculates  true 
AOA  using  both  methods,  and  calculates  the  AOA  position  error 
correction  factor  (from  equation  (2)).  The  program  works  with 
flight  test  data  from  either  USAP  TPS  T-38As  or  RF-4Cs .  No 
data  filtering  is  accomplished. 

Results 

Sample  RF-4C  flight  test  data  from  1  g,  wings  level 
flight  was  processed  by  the  program  AOAOPT.  Test  data  was 
sampled  at  the  highest  rate  possible  for  the  aircraft  DAS,  8 
times  per  second  (a  complete  DAS  description  is  in  Chapter 
IV) .  Figure  3  is  a  plot  of  the  AOA  calculated  using  the  angle 
relationship.  The  calculated  AOA  was  very  sensitive  to  noise 
in  the  altitude  channel  and  was  only  accurate  by  averaging 
over  a  time  span  of  3  to  4  seconds  in  1  g,  wings  level  flight. 


Figure  4  is  a  plot  of  the  AOA  calculated  using  the 

iterative  method.  The  resulting  AOA  is  less  subject  to  noise, 

as  the  pitch  rate  and  normal  acceleration  channels  were  fairly 

noise  free.  The  data  pitch  rate  and  normal  acceleration 

values  were  also  corrected  for  bias  measured  while  on  the 

ground  (bias  was  measured  from  0  deg/sec  pitch  rate  and  1  g 

normal  acceleration) .  The  major  problem  with  the  iterative 

method  is  its  initial  value,  a  ,  which  must  be  calculated 

o 

beforehand.  The  program  AOAOPT  uses  as  aQ  the  AOA  calculated 
from  the  angle  relationship  averaged  over  a  time  interval  of  3 
to  4  seconds  of  1  g  flight.  However,  the  AOA  values  are  still 
corrupted  by  noise  and  are  only  as  good  as  the  resolution  of 
the  DAS.  Some  type  of  filtering  is  needed  to  optimize  the 


Figure  4.  True  AOA  Calculated  Using  the  Iterative  Method  With 
Unfiltered  Data  Compared  to  Measured  AOA  From  the 


III.  State  Estimation 


State  Equations 

In  order  to  use  a  digital  computer  to  filter  the  flight 
test  data  and  compute  an  "optimal"  true  AOA,  the  system 
dynamics  need  to  be  modelled  (6:174).  One  way  to  model  the 
system  is  with  linear  differential  equations  of  the  form: 

x(t)  =  F(t)  x(t)  +  B ( t )  u(t)  (17) 

z(t)  =  H ( t )  x ( t )  (18) 


where  jc  ( t )  is  the  state  vector,  u(t)  is  the  control  input,  and 
z(t)  is  the  measurement  history.  One  differential  equation 
for  angle-of-attack  comes  from  equation  (15): 


a 


q 


(15) 


Another  equation  to  use  in  the  state  equations  is  the  pitch 
rate  equation  valid  for  wings  level  flight: 


8 


=  q 


(19) 


Since  a  and  6  are  not  directly  related,  they  become  functions 
of  the  inputs  q  and  a The  only  useful  parameter  to  measure 
is  pitch  angle,  8,  since  the  measured  AOA  has  an  undetermined 
position  error.  Combining  equations  (15)  and  (19)  together 
gives  the  state  equations: 


“ 
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(20) 
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The  measurement  equation  is: 


W.V T*VvrJV  *7  *  ■ 


[•-]■ co  13  [!] 

The  above  equations  define  the  matrices  F(t),  B(t),  and  H(t). 
This  is  the  math  model  to  compute  a  true  AOA,  but  nothing  in 
the  model  filters  the  noise  in  the  data. 

Kalman  Filtering 

A  Kalman  filter  provides  the  best  method  to  "optimize" 
the  true  AOA  from  available  flight  test  data.  A  Kalman  filter 
will  combine  the  pitch  angle  measurements,  plus  prior 
knowledge  about  the  system  and  measuring  devices,  to  produce 
an  estimate  of  true  AOA  in  such  a  manner  that  the  error  in 
true  AOA  is  minimized  statistically  (7:5).  The  filter  uses 
the  state  equations  plus  a  statistical  description  of  the 
system  noises,  measurement  noises,  and  uncertainty  in  the 
dynamics  model  (7:4).  The  Kalman  filter  assumes  that  the 
system  can  be  described  by  a  linear  model,  and  that  system  and 
measurement  noises  are  white  and  Gaussian  (7:7). 

The  original  system  model,  equations  (17)  and  (18),  is 
augmented  by  (7:146): 

x(t)  =  F ( t )  x(t)  +  B ( t )  u(t)  +  G(t)  w { t )  (22) 

z{t)  =  H ( t )  x(t)  +  v(t)  (23) 

where  the  system  is  now  driven  by  the  input  vector  ij(t)  and 
noise  vectors  w(t)  and  v(t) . 


■» . « m 

V'.'-V 
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The  vector  w(t)  models  the  system  noise  as  white  and 
Gaussian  with  mean  of  zero  and  strength  q,  described  as 
(7:154-155) : 


E[w( t ) ]  =  0  (24) 

E[w(t)  wT(t  +  x)]  =  q  6 ( t )  (25) 

where  q  is  a  measure  of  the  uncertainty  in  the  input  vector 
ia(t).  The  noise  in  the  values  from  the  aircraft  DAS  is 
assumed  to  be  white  since  it  is  random  and  uncorrelated. 

The  vector  v(t)  models  the  measurement  noise  as  white 
and  Gaussian  with  mean  zero  and  strength  r,  described  as 
(7:174) : 


E[v(t)]  =  0  (26) 

E[ v ( t )  vT(t)]  =  R  (27) 

where  R  is  a  measure  of  the  uncertainty  in  the  measurement 
vector  z( t ) . 

In  order  for  the  Kalman  filter  to  propogate  the  system, 

A 

the  estimated  state  vector  (denoted  by  )  must  be  given  an 
initial  condition,  ,  where  (6:233): 

E[x(t  )]  =  x(t  )  (28) 

—  o  —  o 

E[  Cx(tQ)  -  x(to)][x(to)  _  x(to)]T]  =  P(tQ)  (29) 

The  equations  to  propogate  and  update  the  optimal 


estimate  using  Kalman  filtering  are  fully  derived  in 


the  matrix  F(t)  is  a  zero  matrix  (equation  (20)),  the 


propogation  equations  from  a  measurement  at  time  t^_^  to  time 


t .  become : 

i 


x(t.  )  =  xftj^  )  +  At  B(tp  p) 


P(t.  )  =  P(t._1+)  +  At  G(ti)  Q  GT(ti) 


where  denotes  prior  to  the  update  and  denotes  after  the 


update.  The  matrix  J3( t  )  is  defined  in  equation  (20).  Matrix 


G(tp  is  set  equal  to  B(tp  so  the  noise  in  the  input  vector, 


jj(t.),  is  modeled  by  the  values  in  matrix  Q.  Q  becomes  a  2x2 


matrix  which  contains  the  uncertainties  in  the  inputs  q  and 


a^.  These  two  inputs  are  assumed  independent,  therefore  Q 


becomes  a  diagonal  matrix  of  the  form: 


where  the  diagonal  values  are  constant  with  time. 


The  update  equations  at  measurement  time  t^  are: 


K(tp  =  P(tj  )  HT(tp  [  H(tp  P(tj  )  HT(tp  +  R  ]  1 


x(t.  )  =  x(t.  )  +  K(tp  [  Z(t.)  -  Http  x(ti  )  ] 


P(ti  )  =  P(ti  )  -  K(ti)  H ( t i )  P(ti  ) 


The  matrix  H(tp  is  defined  in  equation  (21).  K(tp  is  the 
gain  matrix  which  specifies  how  much  the  measurement,  _Z(tp, 


is  weighted  in  the  update.  R  is  a  lxl  value  modelling  the 


*  »  *  •  f  K  ■  *  *  .  -  . 


.v  .  .s  . 


18 


R  is 


noise  in  the  measurement,  _Z(t.  ),  which  is  pitch  angle, 
in  the  form: 


R 


(36) 


and  is  also  constant  in  time. 

The  Kalman  filter  is  ready  to  be  put  into  a  digital 
computer  routine.  Equations  (30),  (31),  (33),  (34),  and  (35) 

will  propogate  and  update  the  system  over  time.  The  matrix 
_x(t_+)  contains  the  "optimal"  values  for  a  and  e,  of  which  a 
is  the  "optimal"  true  AOA  desired.  Before  the  routine  can  be 
implemented,  the  Kalman  filter  must  be  "tuned"  to  determine 
the  values  for  P(t  ),  Q,  and  R. 


Filter  Tuni ng 

The  objective  of  filter  tuning  is  to  achieve  the  best 
possible  estimation  performance  from  a  filter  that  is  totally 
specified  except  for  _P(tQ),  Q,  and  R.  The  covariance  values 
in  those  matrices  account  for  the  actual  noises  and 
disturbances  in  the  system  and  determine  how  adequately  the 
model  represents  the  real  world  system.  The  P(tQ)  matrix 
determines  the  initial  performance  of  the  filter,  and  the  Q 
and  R  matrices  determine  the  long  term  performance  (7:337). 

The  method  of  filter  tuning  used  here  is  "covariance 
analysis"  (7:337-339).  The  filter  program  is  run  with  some 
assumed  covariance  values  in  the  three  matrices,  P(tQ),  Q* 
and  R.  The  "true”  root  mean  squared  ( rms )  error,  which  is  the 
error  at  each  update  between  the  filter's  estimate  and  the 


actual  measurement,  is  plotted  over  time.  The  time  history  of 
the  computed  rms  error,  or  what  the  filter  calculates  as  its 
error,  is  plotted  with  the  true  error. 


Figure  5 


M 


C  nmpntctl  riih  error 
True  "  rms  error 


Measurement  Weighed 
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An  Example  of  Kalman  Filter  Tuning  Through 
Covariance  Analysis  (7:338) 


Figure  5  is  an  example  of  what  these  plots  show.  Plot 
5(a)  shows  a  filter  that  has  a  low  computed  error  and  weighs 


the  measurement  too  little.  Plot  5(b)  shows  a  filter  that  has 

too  high  a  computed  error  and  weighs  the  measurements  too 

much.  Plot  5(c)  shows  a  filter  that  is  just  right  -  its 

computed  error  and  true  error  are  equal  (7:339).  For  the 

Kalman  filter  used  in  this  project,  the  values  of  P(t  ),  Q, 

—  o  — 

and  R  were  varied  until  the  computed  rms  errors  and  true  rms 
errors  were  about  equal.  The  Kalman  filter  was  complete  and 
ready  to  filter  test  data. 

Computer  Program  KALOPT 

A  FORTRAN  computer  program  was  designed  to  calculate  an 

"optimal"  true  AOA  using  the  Kalman  filter  equations.  This 

program  is  called  KALOPT  and  is  shown  in  Appendix  C.  The 

program  reads  the  required  flight  test  parameters  from  the 

DAS,  makes  necessary  pitot-static  corrections,  and  calculates 

an  "optimal"  true  AOA  each  iteration.  The  program  works  with 

flight  test  data  from  either  USAF  TPS  T-38As  or  RF-4Cs . 

The  Aydin-Vector  DAS  does  not  read  normal  acceleration, 

a  ,  but  instead  reads  normal  load  factor,  n  .  The  sign  of  n 
z  z  3  z 

is  opposite  from  the  standard  body  axis  system:  positive  n^ 

is  through  the  top  of  the  canopy.  Also,  n z  includes 

acceleration  due  to  gravity.  The  following  equation,  which 

assumes  small  pitch  and  roll  angles,  corrects  n  to  a  : 

z  z 

a,  =  -  ( n  -  1)  32.2  (37) 

z  z 


The  Kalman  filter  needs  to  know  the  initial  values  for 
AOA  and  pitch  angle  to  use  as  x(t  ).  KALOPT  uses  the  first 


DAS  value  for  pitch  angle  as  eQ.  However,  an  initial  AOA 
needs  to  be  calculated  since  the  DAS  AOA  values  have  the 
yet-to-be-determined  position  error.  The  value  of  aQ  is 
calculated  using  the  angle  method  (equations  (1)  and  (5))  used 
in  the  program  AOAOPT.  It  averages  the  unfiltered  true  AOAs 
over  a  3  to  4  second  period  in  1  g  flight  to  calculate  aQ . 

KALOPT  calculated  the  computed  and  true  rms  errors  after 
each  iteration.  The  values  were  varied  from  0.001  to  1.0 
during  the  covariance  analysis.  Changing  the  values  of  P(tQ) 
changed  the  initial  rms  values,  but  had  little  effect  on  the 
overall  results.  When  the  Q  values  were  increased,  the 
measurement  was  weighted  more;  the  true  rms  error  was  less 
than  the  computed  rms  error.  Increasing  the  R  value  caused 
the  measurement  to  be  weighted  less;  the  true  rms  error  was 
greater  than  the  computed  rms  error.  These  results  agreed 
with  the  theory  behind  filter  tuning .  Based  on  the  covariance 
analysis  conducted  using  T-38A  data,  the  following  values 
caused  the  true  and  computed  errors  to  be  equal; 

P(t  )  .  r  °-io°  °  i 

-  o'  0  0.030 

q  -  [  °-025  0  1 

-  0  0.002 

R  =  [0. 300 

KALOPT  reads  in  the  covariance  values  from  a  separate  file,  so 
they  can  be  easily  changed  without  changing  the  program. 


Results 

The  first  attempt  at  Kalman  filtering  included  another 
state  equation  formed  by  combining  equations  (1)  and  (5): 

h  =  (  e  -  a)  U 

The  DAS  altitude  readout  was  used  as  a  measurement  along  with 
pitch  angle.  Unfortunately,  the  noise  of  the  altitude 
transducer  in  the  T-38A  DAS  was  too  erratic  and  could  not  be 
modelled  as  Gaussian.  Altitude  was  not  used  in  KALOPT. 

The  program  KALOPT  processed  the  same  RF-4C  flight  test 
data  that  was  used  in  Chapter  II.  The  test  data  was  sampled  8 
times  per  second.  Figure  6  is  a  plot  of  the  "optimal"  true 
AOA  calculated  by  KALOPT  from  that  data.  The  measured  AOA  and 
the  unfiltered  true  AOA  calculated  by  AOAOPT  are  also  shown. 
The  "optimal"  true  AOA  is  quicker  to  return  to  a  steady  state 
value  than  the  unfiltered  true  AOA.  It  is  impossible  to  tell 
which  AOA  is  more  accurate  as  the  actual  true  AOA  is  unknown. 

A  better  way  to  see  how  the  Kalman  filter  is  working  is 
to  compare  the  measured  pitch  angle  to  the  "optimal"  pitch 
angle  to  see  how  well  it  filters  over  noise  and  resolution 
increments.  Figure  7  is  a  plot  of  measured  pitch  angle  and 
"optimal"  pitch  angle.  The  measured  pitch  angle  only  had  a 
resolution  of  0.7  degrees,  and  after  two  samples  it 
immediately  increased  by  that  amount.  The  "optimal"  pitch 
angle  is  a  fairly  smooth  curve  over  the  time  span,  which  shows 
that  the  Kalman  filter  is  working.  The  next  step  is  to  use 
flight  test  data  to  calculate  the  AOA  position  error. 


The  test  aircraft  was  a  USAF  TPS  T-38A  Talon.  The  T-38A 


is  a  two  place  (tandem)  jet  trainer  which  is  used  extensively 
in  the  TPS  curriculum.  The  aircraft  is  powered  by  two 
J85-GE-5  turbojet  engines  which  give  it  a  maximum  capability  of 
Mach  1.2  in  level  flight  (8:6-6).  Figure  8  is  a  photograph  of 
a  T-3RA  used  at  the  TPS.  A  single  T-38A,  serial  number  (S/N) 
68-8205,  was  used  for  all  data  flights  in  this  project. 
Statistics  on  that  airplane  are  shown  in  Table  I.  An 
important  measurement  is  the  distance  from  the  eg  of  the 
aircraft  to  the  AOA  measuring  vane,  25  feet.  This  length  (x^) 


Figure  8.  USAF  T-38A  Talon 


is  used  to  correct  the  true  AOA  for  pitch  rate  (see  equation 
(3)).  The  eg  of  the  T-38A  only  shifts  0.3%  mean  aerodynamic 
chord  (MAC)  while  consuming  fuel,  which  is  only  0.25  inches. 
Therefore,  the  Ac g  term  from  equation  (3)  can  be  neglected. 


TABLE  I 

USAF  T-38A  Talon  Statistics  (S/N  68-8205)1 


Engines:  Two  J85-GE-5  Turbojets 
Dimensions : 

Length:  46  ft 
CG  to  AOA  Vane:  25  ft 
Wingspan:  25  ft 
MAC :  7  ft 
Height:  12  ft 


Weights : 

Operating  Weight: 

Fuel  (JP-4)  Weight: 
Takeoff  Gross  Weight: 

8.  533 
3,790 
12, 323 

lbs 

lbs 

lbs 

Center  of  Gravity  Movement: 

CG  w/  3,790  lbs  fuel: 

18.4% 

MAC 

CG  w /  400  lbs  fuel : 

18.1% 

MAC 

(8:1-1) 

The  test  T-38A,  S/N  68-8205,  was  modified  for  flight 
testing.  The  most  important  modification  concerning  this 
project  is  a  fully  instrumented  yaw  and  pitch  system  (YAPS) 
noseboom  (9:A.l),  shown  in  Figures  9(a)  and  9(b).  A  complete 
diagram  of  the  YAPS  noseboom  is  in  Appendix  E.  The  YAPS  nose 
boom  has  two  vane-type  sensors,  one  for  AOA  and  one  for 
sideslip  angle.  These  vanes  are  in  front  of  the  fuselage, 


4  in 
0  in 
3  in 
0  in 
11  in 


away  from  the  aerodynamic  influence  of  the  aircraft,  so  the 
AOA  (and  sideslip)  position  errors  should  be  less  than  for 
fuselage  mounted  sensors. 

The  YAPS  noseboom  on  the  T-38A  is  made  out  of  aluminum 
alloy.  It  has  been  structurally  tested  up  to  8.3  g  and  only 
minimal  bending  resulted  ( 9 : D. 43-D . 53 ) .  As  a  result,  bending 
was  ignored  during  this  evaluation.  The  YAPS  noseboom  is 
canted  4  degrees  down  from  the  aircraft  centerline. 

Instrumentation 

An  internal  Aydin-Vect.or  SAU-537  DAS  was  installed  in 
the  aircraft  to  measure  flight  test  parameters.  The  following 
components  of  the  DAS  were  used  in  this  project:  a  vertical 
gyro  installed  in  the  nose  section  to  measure  pitch  and  roll 
angles;  a  three  axis  rate  gyro  installed  in  the  nose  section 
to  measure  pitch,  roll,  and  yaw  rates;  a  three  axis 
accelerometer  installed  in  the  center  fuselage  (at  the  nominal 
eg  location)  to  measure  acceleration  in  the  x,  y,  and  z  axes 
(10:1.1-1.8).  Other  instruments  were  installed  in  the 
aircraft  for  flight  test,  but  they  were  not  used  in  this 
project . 

The  test  aircraft  was  equipped  with  an  internal  Conrac 
ATR-580T70  magnetic  tape  recorder  in  the  aft  cockpit  to  record 
the  data  parameters  (9:A.l).  Forty-eight  data  channels  were 
recorded.  Indicated  airspeed  and  altitude  were  recorded  with 
16  bit  precision,  the  other  parameters  for  this  project  had  8 
bit  precision.  Table  II  is  a  summary  of  the  parameters  used 
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in  this  project  with  their  maximum/minimum  values,  precision 
and  accuracy.  The  parameters  were  recorded  8  times  per 
second . 


Parameter 

Units 

Min 

Value 

Max 

Value 

Resolution 

Accuracy 

Altitude 

feet 

0 

65000 

1.030 

0.103 

Airspeed 

knots 

0 

1250 

0.019 

0.0019 

AOA 

degrees 

-22 

28 

0.202 

*0.101 

Sideslip 

degrees 

-20 

20 

0.164 

0.082 

Pitch 

degrees 

-80 

80 

0.704 

0.704 

Roll 

degrees 

-180 

180 

1.408 

2.816 

Pitch  Rate 

deg/sec 

-20 

20 

0.163 

0.163 

N 

z 

g 

-3 

6 

0.037 

0.0037 

*  Actual  accuracy  +/-  0.5  degrees  due  to  hysteresis 


Before  any  flight  test  was  performed,  all  of  the  DAS 
instruments  were  ground  calibrated  and  their  calibration  files 
updated.  All  important  instruments  were  found  to  be  working 
correctly  except  for  the  AOA  transducer.  It  had  a  large 
hysteresis  problem  due  to  wear  on  its  internal  gearing.  This 
hysteresis  is  shown  in  Figure  10.  There  is  a  +/-  0.5  degree 
error  in  true  AOA  depending  on  whether  the  vane  is  moving  up 
or  down.  The  AOA  transducer  was  designed  in  the  1960s,  and  no 
replacement,  parts  are  available.  The  worn  gears  could  not  be 
fixed  or  replaced.  This  hysteresis  will  have  a  large  effect 
on  the  flight  test  data. 


Flight  Test  Method 

The  purpose  of  the  flight  testing  was  to  gather  data  to 
calculate  true  AOA  and  to  see  what  factors  affected  the  AOA 
position  error.  AOA  position  error  is  primarily  a  function  of 
AOA,  therefore  the  testing  covered  large  AOA  changes.  Other 
possible  factors  that  were  considered  in  designing  the 
maneuvers  were  Mach  number,  Reynold's  number,  and  sideslip 
angle.  All  testing  was  conducted  wings  level  due  to  the 
assumptions  used  in  the  AOA  equations  (see  Chapter  II). 

Since  Reynold's  number  was  a  possible  factor,  testing 
was  performed  at  different  altitudes.  Due  to  the  altitude 
restrictions  of  available  supersonic  airspace,  25,000  feet  and 
15,000  feet  were  chosen  for  the  testing.  In  order  to  see  the 
effects  of  sideslip  on  AOA,  the  first  maneuver  to  be  performed 
was  a  wings  level,  slowly  varying  sideslip  using  maximum 
rudder  deflection  in  both  directions.  Thrust  was  varied  to 
maintain  Mach  number  constant.  This  maneuver  was  performed  at 
different  Mach  numbers.  Actual  data  points  are  shown  in 
Appendix  D  (11:3). 

After  the  sideslip  maneuver,  a  roller  coaster  maneuver 
was  performed  to  vary  AOA  as  much  as  possible.  From  a  1  g 
trim  condition,  the  aircraft,  nose  was  pulled  up  slightly,  then 
pushed  forward  to  the  minimum  load  factor  specified  for  that 
data  point.  An  onset  rate  of  3  seconds  per  g  minimum  was 
desired  throughout  the  maneuver.  At  the  minimum  load  factor, 
the  aircraft  nose  was  pulled  back  to  the  maximum  load  factor 


specified  for  that  data  point.  The  aircraft  nose  was  then 
pushed  forward  to  regain  1  g  level  flight.  All  data  points 
and  load  factor  limits  are  in  Appendix  D  (11:3).  Thrust  was 
varied  to  maintain  constant  Mach  number  during  the  maneuver. 

All  testing  was  performed  in  the  cruise  configuration 
(gear  and  flaps  up)  with  no  external  stores.  All  T-38A  FI ight 
Manual  (8)  limitations  were  complied  with.  Additional 
restrictions  in  the  T-38A  AOA  Position  Error  Test  Plan  (11) 
were  followed. 

Test  Results 

Three  T-38A  test  flights  were  flown  at  the  USAF  Flight 
Test  Center,  Edwards  AFB ,  California.  A  summary  of  these 
flights  is  shown  in  Appendix  D.  No  data  were  gathered  on  one 
flight,  due  to  bad  weather.  The  same  aircraft,  68-8205,  was 
used  on  all  three  flights  due  to  scheduling  availability. 
Future  test  programs  using  this  method  should  fly  different 
tail  numbers  to  prevent  bias  from  one  aircraft's  own 
peculiarities . 

The  wings  level  sideslip  maneuver  was  performed  at  all 
data  points.  The  maximum  sideslip  angle  generated  was  +/-  3.8 
degrees  at  25,000  feet  pressure  altitude  (Hc),  Mach  0.45.  No 
change  in  AOA  was  found  at  this  point  or  any  of  the  others. 

In  wind  tunnel  testing  performed  on  a  Conrac  NB1U,  a  noseboom 
similar  to  the  T-38A  YAPS  noseboom,  no  AOA  position  error 
change  was  discovered  until  five  degrees  of  sideslip 
(3:38-39) . 
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The  roller  coaster  maneuver  was  performed  at  all  data 
points,  and  repeated  at  the  25,000  feet  Hc  points.  The  data 
was  reduced  using  the  FORTRAN  program  KALOPT  (see  Chapter 
III).  Measured  AOA  was  plotted  against  the  "optimal"  true  AOA 
(henceforth  referred  to  as  true  AOA)  at  each  point  tested. 
These  plots  are  Figures  12  -  29,  Appendix  A.  The  data  points 
plot  out  fairly  linear,  which  shows  that  the  equations  and 
Kalman  filtering  worked.  Most  of  the  test  points  flown  up  to 
6  g  were  terminated  at  that  point  due  to  the  Mach  number 
decreasing  outside  tolerances  (+/-  0.02  Mach  desired).  Also, 
at  many  of  the  high  g  points  the  data  trace  becomes  erratic. 
This  was  due  to  aerodynamic  buffet.  Future  test  maneuvers  for 
this  method  do  not  have  to  go  to  such  high  g  limits,  as  the 
data  collected  at  lower  g  limits  is  satisfactory. 

As  expected,  the  hysteresis  error  due  to  mechanical  lag 
in  the  AOA  gears  was  evident  in  the  results.  All  eighteen 
plots  show  two  lines  of  data,  depending  on  whether  the  AOA 
vane  was  moving  down  or  up.  The  error  between  the  two  lines 
ranges  from  +/-  0.5  to  +/-  0.8  degrees,  similar  to  the 
hysteresis  error  in  Figure  10.  In  order  to  average  the  error, 
a  straight  line  was  drawn  down  the  middle  of  the  two  lines. 

The  slope  of  this  line  is  the  AOA  position  error  correction 
factor,  K ^  (from  equation  (2)).  The  x-axis  intercept,  a°T» 
was  also  determined  from  these  plots.  These  values  are 
summarized  in  Table  III.  The  Reynold's  numbers  were 
calculated  using  MAC  (7  feet)  as  the  constant  length. 


TABLE  II 


Summary  of  Flight  Test  Results 


Date  of  Flight 

A1 1 it  ude 

Mach 

Reynolds  # 

K 

a 

p 

0 

|-3 

(  feet ) 

(x  107) 

(deg ) 

25,000 

0.84 

1.99 

1 . 26 

1 .25 

25,000 

0.94 

2.23 

1.38 

0.60 

25,000 

0.98 

2. 30 

1 . 32 

2.00 

25,000 

1.07 

2.51 

1.48 

2.30 

25,000 

0.62 

1 . 46 

1 .  16 

1.20 

25,000 

0.44 

1.03 

1 .  16 

1.20 

■lilliM 

15,000 

0.83 

2.69 

1  .  10 

m 

15,000 

0.94 

3.03 

1 .14 

15,000 

0.96 

3.  10 

1. 33 

mu 

15,000 

1.07 

3.44 

1.45 

0.90 

15,000 

0.64 

2.08 

1.30 

0.80 

HHH 

15,000 

0.42 

1.35 

1.20 

1.90 

25,000 

0.81 

1.90 

1 .  10 

1 .  10 

25,000 

0.92 

2.16 

1.20 

2.00 

25,000 

0.96 

2.24 

1.30 

*-1.9 

25,000 

1.06 

2.47 

1.36 

2.20 

25,000 

0.65 

1.53 

1 .16 

0.70 

WwB 

25,000 

0.44 

1.03 

1 . 20 

0.50 

*  Exceeds  2  standard  deviations  from  mean 


The  values  for  K  are  plotted  versus  Mach  number  in 

a 

Figure  11.  A  curve  was  drawn  through  the  points  and  shows  a 

large  increase  in  K  as  Mach  number  increases  above  0.8.  No 

a 

apparent  Reynold's  number  effects  in  K  are  evident  in 

a 

comparing  the  15,000  feet  points  to  the  25,000  feet  points. 

The  K  values  from  the  Conrac  NBIU  wind  tunnel  testing  are 
a 

also  plotted  in  Figure  11  (3:38).  Although  the  shapes  of  the 

curves  are  similar,  the  values  for  K  are  different. 
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Vought  YAPS  Noseboom 


both  nosebooras  are  so  similar  in  size  and  shape,  the 
difference  is  mainly  due  to  the  lack  of  fuselage  and  wing 
effects  on  the  wind  tunnel  results. 

Equation  (2),  the  equation  to  calculate  true  AOA,  is: 


m 


(  a  + 


K 


U 


\  a ) 


(2 


U 


For  T-38A  68-8205,  (x  +  Ac g)  is  assumed  a  constant  25  feet. 

a 

v  is  6  inches  (see  Appendix  E) .  The  values  for  K  are  shown 
•'a  a 

in  Figure  11  as  a  function  of  Mach  number.  However,  equation 
(2)  assumes  that  no  AOA  position  error  exists  at  zero  degrees 
AOA.  According  to  Figures  12  to  29,  this  is  not  true  for  the 
T-38A.  Some  bias  exists,  which  is  the  x-axis  intercept,  aoT. 
Adding  this  bias  to  equation  (2),  and  neglecting  pitch  and 
roll  rate,  gives: 


a 


c 


K 

a 
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The  values  for  aoT  from  Figures  12  to  29  are  shown  in  Table 
III.  The  values  vary  randomly,  and  do  not  seem  to  be 
functions  of  Mach  number.  The  average  of  all  18  values  is 
1.26  degrees,  with  a  standard  deviation  of  1.06  degrees.  One 
value,  for  25,000  feet  Hc  and  Mach  0.96,  is  -1.9  degrees, 
which  exceeds  two  standard  deviations  from  the  mean. 


Neglecting  that  point  as  erroneous,  the  average  of  the 


remaining  values  is  1.45  degrees,  with  a  standard  deviation  of 
0.73  degrees.  The  equation  to  solve  for  the  true  AOA  of  the 
T-38A  with  a  Vought  YAPS  noseboom  is: 


1.45 


m 


The  pitch  and  roll  rate  terms  from  equation  (2)  should  be 
included  when  applicable. 


(39) 


V.  Conclusions  and  Recommendat ions 

All  project  objectives  were  met.  Conclusions  and 
recommendations  follow  in  order  of  importance: 

The  Kalman  filter  program  KALOPT  calculated  "optimal" 
true  angle-of-attack  ( AOA)  values  for  a  T-38A  Talon  using 
equations  of  motion  for  wings  level  flight  and  pitch  angle 
measurements.  From  these  true  AOA  values,  the  AOA  position 
error  correction  factors  were  determined  and  were  found  to  be 
functions  of  Mach  number.  Only  standard  T-38A  flight  test 
instrumentation  was  used  -  no  flight  path  accelerometers  were 
needed.  This  method  proved  to  be  accurate  in  gathering  data 
with  minimal  instrumentation  over  a  large  range  of  AOAs . 

1.  THIS  STATE  ESTIMATION/KALMAN  FILTERING  METHOD  OF 
CALCULATING  AOA  POSITION  ERROR  SHOULD  BE  CONSIDERED  IN 
FUTURE  AOA  ERROR  TESTING. 

The  Kalman  filter  needed  an  initial  AOA  to  start 
propogating  an  "optimal"  true  AOA.  This  initial  AOA  was 
calculated  from  1  g  wings  level  flight  prior  to  the  roller 
coaster  flight  test  maneuver.  Also,  aerodynamic  buffet  at 
high  load  factors  caused  some  data  scatter.  The  data  gathered 
prior  to  the  buffet  were  enough  to  calculate  the  AOA  position 
error . 

2.  FUTURE  MANEUVERS  TO  GATHER  DATA  FOR  THIS  METHOD  SHOULD 
START  FROM  A  1  G  TRIM  SHOT  FOR  3  TO  4  SECONDS.  THE 


MANEUVERS  SHOULD  TERMINATE  PRIOR  TO  AERODYNAMIC  BUFFET 


The  T-38A  flight  test  data  were  not  completely  accurate 
due  to  hysteresis  errors  in  the  AOA  transducer.  The  accuracy 
of  the  method  would  be  better  determined  using  an  aircraft 
with  no  AOA  hysteresis  error.  Also,  testing  with  an  aircraft 
with  a  fuselage-mounted  AOA  sensor  would  show  larger  AOA 
position  errors  and  would  further  validate  the  method. 

3.  FURTHER  TESTING  SHOULD  BE  CONDUCTED  USING  AN  RF-4C  OR 
OTHER  SUITABLE  AIRCRAFT  THAT  HAS  NO  AOA  HYSTERESIS  ERROR 
AND  HAS  A  FUSELAGE-MOUNTED  AOA  SENSOR. 

4.  THE  USAF  TEST  PILOT  SCHOOL  NEEDS  TO  INSTALL  NEW  AOA 
TRANSDUCERS  IN  THEIR  T-38A  AIRCRAFT  TO  ELIMINATE  THE  AOA 


HYSTERESIS  ERRORS. 
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Figure  17.  T-38A  Angle-of-Attack  Position  Error  for  a  Vought 

YAPS  Noseboom  AT  25,000  feet  H  ,  Mach  0.44 
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Figure  23.  T-38A  Angle-of-At.tack  Position  Error  for  a  Vought 

YAPS  Noseboom  at  15,000  feet  H  ,  Mach  0.42 
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PROGRAM  AOAOPT 

C 

C  PROGRAM  AOAOPT  IS  A  FORTRAN  PROGRAM  DESIGNED  TO  CALCULATE 
C  ANGLE-OF-ATTACK  FROM  PARAMETERS  RECORDED  INFLIGHT  ON  AN 
C  AYDIN-UECTOR  DATA  ACQUISITION  SYSTEM.  AOAOPT  USES  TIME  (TI), 

C  IAS  (AS).  INDICATED  ALT  (ALT).  MEASURED  AOA  (AOA) >  PITCH 
C  ANGLE  (TH),  PITCH  RATE  (Q),  NORMAL  ACCELERATION  (ZN).  SIDESLIP 
C  ANGLE  (AOSS).  AND  BANK  ANGLE  (BA)  OBTAINED  FROM  DYNAMICS  EUP 
C  AND  EUS  FILES.  AOAOPT  CORRECTS  FOR  PITOT-STATIC  ERRORS  TO 

C  DETERMINE  PRESSURE  ALTITUDE  (HC).  MACH  NUMBER  (AMC).  AND 
C  TRUE  AIRSPEED  (UTAS).  THE  PROGRAM  CALCULATES  TRUE  AOA  BY  TWO 
C  METHODS i  ANGLE  RELATIONSHIP  WHERE  AOA  -  PITCH  -  FLIGHT  PATH. 

C  AND  ITERATIUE  METHOD  WHERE  AOA-DOT  -  PITCH  RATE  -  NORMAL  G. 

C  PROGRAM  CAN  BE  USED  WITH  T-38A  OR  RF-4C  DATA, 

e 

BYTE  PNAMEC 12)  .ENAME(  12)  >  DOF(?)  ■  ATYPE(S) .  ATAILO) 

DIMENSION  TI  (300)  .AS  (300) .  ALT  (300) .  AOA  (300) .  TH(300)  »Q(300)  , 

1ZN ( 300 ) . AOSS (300 ) . BA ( 300 ) ,HC ( 300 ) , AMC (300 ) , UTAS ( 300 ) . A0A2( 300 ) . 
2A0A3 ( 300 ), HDOT ( 300 ). FPANG ( 300 ), AOA 1 ( 300 ). ERROR ( 300 )• ADOT (300 ) . 
3AD0TM  ( 300 ),  ERROR 1  ( 300 ).  ERR0R2  ( 300 ) 

BYTE  FILEl(lS)  .OFILEK17) 

BYTE  FILE2(15) . 0FILE2( 17) 

BYTE  FILE3C IS) ,0FILE3( 17) 

BYTE  FILE4(15) «0FILE4( 17) 

BYTE  FILE6(15),0FILE6(1?) 

DATA  FILEV’E’ , ’S' . *  »’ .0*0, ’ . ' , ’E* , ’U’ , 'PV 
DATA  FILE2/’ E’ , ’S' , ’ 1 ’ ,8*0, ’  .  ’  ,  ’ E’ , ’ U’  .  ' S  V 
DATA  FILE3" E’ , ’S’ , ’ «’ ,0*0, ’ . ’ , ’D’ . ’A’ , ’TV 
DATA  FILE4/’E’ , ’S’ . ’ ■’ ,8*0, ’ . ’ , ’D’ , ’A’ , ’T’/ 

DATA  FILE6."E’ .  ’S’ .  ’  '  *  ,8*0,  ’ .  ’ ,  'D' ,  'A' ,  'TV 
C 

WRITE  (5,100) 

100  FORMAT ( ’  ENTER  NAME  OF  DYNAMICS  FILE'  NNNNN . EUP ’ , / ) 

ACCEPT  110,  (FILEKI).  1-4,11) 

110  FORMAT  (BA1) 

OPEN(UNIT-l,NAME-FILEl,TYPE=’UNKNOWN’ ) 

WRITE  (5,120) 

120  FORMAT ( ’  ENTER  NAME  OF  DYNAMICS  FILE'  NNNNN. EUS’ , /) 

ACCEPT  110, (FILE2( I ), 1-4,11) 

OPEN(UNIT-2,NAME-FILE2, TYPE- ’UNKNOWN’ ) 

WRITE(S,13v 

130  FORMAT ( ’  ENTER  NAME  FOR  DATA  FILE  1'  NNNNNN . DAT ’ , / ) 

ACCEPT  110,  (FILE3H).  1-4,11  ) 

OPEN ( UN I T -3 , NAME - F I LE3 , TYPE - ’ UNKNOWN ’ ) 

WRITE(S,  14f 

140  FORMAT ( ’  ENTER  NAME  FOR  DATA  FILE  2=  NNNNNN . DAT ’ > / ) 

ACCEPT  110, (FILE4( I). 1-4,11 ) 

OPEN  (  UN  I T  -  4 ,  NAME -FIL.E4,  TYPE  -  ’  UNKNOWN  ’  ) 

WRITE(S,  15V 

150  FORMAT (’  ENTER  NAME  FOR  DATA  FILE  3 '  NNNNNN. DAT’ , /) 

ACCEPT  110, (FILE6( I ), 1-4,11 ) 
0PEN(UNIT-6,NAME-FILE6,TYPE=’UNKN0WN’ ) 

WRITE(5, IGA' 

160  FORMAT  (’  ENTER  NUMBER  OF  DATA  POINTS  DESIRED  '  ,/) 

READ(5. 170)  N 
170  FORMAT  (IS) 


WRITE(5> 190) 

FORMAT  (’  ENTER  FIRST  LINE  OF  DATA  DESIRED  ’  ./) 

READCSf 190)  M 
FORMAT  (15) 

WRITE(5.192) 

FORMAT  (’  ENTER  NUMBER  OF  INITIAL  CONSTANT  AOA  LINES  ’  ./) 

READ(5.19S)  NA 

FORMATCIS) 

WR!TE(S,200> 

FORMAT ( *  ENTER  CORRECTIONS  FOR  PITCH  RATE  AND  NZ»  XXX.YYY’ ./. 

1’  CORR  ARE  FROM  GROUND  BLOCK,  ABOUE/BELOW  0  FOR  Q’,/, 

2’  ABOUE/BELOW  1  FOR  NZ’ ,/) 

READ (5,210) QCORR , ZNCORR 
FORMAT (2F10. 3) 

THIS  PORTION  READS  DATA  FROM  THE  DYNAMICS  EUP  AND  EUS  FILES 
AND  FORMATS  THREE  DATA  FILES  (XXXX.DAT)  WHERE  THE  PROGRAM 
RESULTS  WILL  BE  SENT.  DATA  FILE  #1  RECORDS  THE  RAW  DATA  FROM 
THE  EUP  AND  EUS  FILES.  DATA  FILE  42  CONTAINS  AOA  COMPUTED 
BY  THE  ANGLE  METHOD.  DATA  FILE  43  CONTAINS  AOA  COMPUTED  BY 
THE  ITERATIVE  METHOD. 

READ  (1,220 )PNAME,ENAME, DOF 
FORMAT  (8X.12A1,15X,12A1.21X,?A1) 

WRITE  (3,230)PNAME«ENAME,DOF 
WRITE  (4, 230)PNAME, ENAME, DOF 
WRITE  (6,230)PNAME,ENAME, DOF 

FORMAT  (’  PILOT ' ’ , IX, 12A1 , SX, *  ENGINEER*  * , IX, 12A1.5X. 
l’DATE  OF  FLIGHT'’ , IX, 7A1) 

READ  (1,240)ATYPE> ATAIL 
FORMAT  (llX.SAl, 13X.3A1) 

WRITE  (3,2S0)ATYPE,ATAIL 
WRITE  (4,2S0)ATYPE> ATAIL 
WRITE  (6,250)ATYPE, ATAIL 

FORMAT ( ’  A/C  TYPE' ’, IX, 5A1,5X, 'TAIL  4«’,1X,3A1) 

WRITE  (3,260) 

FORMAT  ( // , TX ,  ’LINE' ,6X, ’TIME’ ,2X, ’AIRSPEED’ ,2X, ’ALTITUDE’ ,7X. 
l’A0A’,5X, ’PITCH’ ,2X, ’PITCH  RT’ ,8X, ’NZ’ ,6X, ’AOSS’ ,2X, ’ROLL  ANG’ ) 
WRITE  (4,270) 

FORMAT  (//. 7X,’ LINE’, 6X, ’ TIME’, IX, ’PRESS  ALT’ ,6X, ’MACH’ ,3X, ’TRUE 
1  AS’.TX. ’R0C’,2X, ’FLT  PATH’ ,5X, ’AOA-l’ ,SX, ’AOA-M’ ,5X, ’ERROR’ ) 
WRITE  (6,280) 

FORMAT  (//,7X, ’LINE’ ,6X, ’TIME’ ,4X, ’A-DOT1’ ,4X, ’ A-DOTM’ ,SX, ’AOA-2 
15X, ’AOA-M’ ,5X, ’ERROR’ ,5X, ’AOA-3’ ,5X, ’ AOA-T’ ,5X, ’ERROR’ ,3X, 
2'K-ALPHA’ ) 

WRITE  (3,290) 

FORMAT  ( 16X, ’ (SEC)’ ,3X. ’ (KNOTS)’ ,4X, ’ (FEET)’  ,5X,  ’ (DEG)’ , 

15X, ’ (DEG) ’ , IX, ’ (DEG/SEC) ' , 7X, ’ (G) ’ ,SX, ’ (DEG) ’ ,5X, ’ (DEG) ’ ,/) 

WRITE  (4,300) 

FORMAT  ( 16X, ’ (SEC)’ ,4X, ’ (FEET)’ , 12X, ’ (FT/SEC) ’ ,2X, 

1* (FT/SEC)’ ,5X, ’ (DEG)’ ,5X, ’ (DEG)’ ,5X, ’ (DEG)’ ,5X, ’ (DEG)’ ./) 
WRITE(6,310) 

FORMAT  ( 16X,  ’  (SEC)  ’ ,  IX,  ’  (DEG/SEC)  MX,  ’  (DEG/SEC)  ’  ,5X,  ’  (DEG)  ’ , 
1SX, ’ (DEG) ’ ,5X, ’ (DEG) ’ ,5X. ' (DEG) ’ ,SX, ’ (DEG) ’ ,5X, ’ (DEG) ’ ,/) 

READ  (1,320) 

FORMAT  (/,/,/) 


REF®  (2.330) 

330  FORMAT  (/,/. /,/,/) 

DO  350  J-l.M 
READ(1,340) 

READ(2»340) 

340  FORMAT  (  ) 

350  CONTINUE 

DO  390  1*1.  N 

REF®(1, 3601TI  ( I ) ,  AS(I ) .  FlLTC  I)  ,AOAd ) .  TH<  I ) .  QC I ) .  ZN<  I ) 

360  FORMAT  (31X,F10.3,F10.3,F10.3,F10.3,F10.3,F10.3,20X,F10.3) 

READ(2, 370) AOSS(  I ) ,  BAC I ) 

370  FORMAT  (31X.F10.3. 10X.F10.3) 

C 

C  THE  NEXT  TWO  STEPS  CORRECT  PITCH  RATE  AND  NORMAL  G  FOR  BIAS 
C  FOUND  ON  THE  GROUND  (BIAS  COMPUTED  FROM  0  DEG/SEC  FOR  PITCH 
C  RATE  AND  1  G  FOR  NORMAL  G) 

C 

Q(I)»Q(I)-QCORR 
ZN( I )»ZN( I )-ZNCORR 

WRITE(3, 380) I . TI (I ) . AS( I ) . ALT(I ) , AOA( I ) . TH( I ) .Q( I ) ,ZN( I ) , AOSS( I ) , 
1BA(I) 

380  FORMATdX.  U0.9F10.3) 

390  CONTINUE 

C 

C  THIS  PORTION  PERFORMS  THE  PITOT-STATIC  CORRECTIONS  TO  COMPUTE 
C  PRESSURE  ALT,  MACH  NUMBER,  AND  TRUE  AIRSPEED.  IT  USES  THE 
C  STANDARD  PITOT-STATIC  EQUATIONS  AND  ERROR  COEFFICIENTS  FROM  THE 
C  USAF  TEST  PILOT  SCHOOL  FILES. 

C 

DO  710  I-l.N 

IF(ALT(I)  .GT.  36089) GO  TO  500 
DELTA* (1 - ( 6 . 87559E-6*ALT ( I ) ) ) **5.2561 
THETA* 1- ( 6 . 87S59E-6*ALT ( I ) ) 

GO  TO  510 

500  DELTA- . 22337*EXP ( -4 . 80635E-S* ( ALT ( I ) -36089 ) ) 

THETA*. 751874 

510  IF(AS<I)  .GT.  661 . 48)G0  TO  520 

QCICPA*(d+(.2*((AS(I)/661.48)**2)))**3.5)-l 
GO  TO  530 

520  QCICPA-((166.922*((AS(I)/661.48)**7))/(((7*((AS(I)/661.48)**2)) 

1-1)**2.S))-1 
530  QC I CPS - QC I CPA/DELTA 

AMIC-SQRT(5X( ( (QCICPS+1 )**( . 2857143) )-l ) ) 

IF(ATYPEd)  .EQ.  ’T’)GOTO600 
C 

DHPC- ( -907+ (10270*AM I C )  +  ( -44495* < AM I CX*2 ) ) +  < 93931* ( AM I  DM3 ) ) 

1 + ( -95540* ( AM I C**4 ) ) + ( 37208* ( AM I C**5 ) ) )*THETA 
IF(AMIC  .GT.  1)DHPC=0 
GO  TO  700 

600  IF(AMIC  .LT.  . 955)G0  TO  620 
IF(AMIC  .LT.  . 967) GO  TO  640 
IF(AMIC  .LT.  1 . 02S)GO  TO  660 
C0— 46325 
Cl -82583 
C2— 36667 
GO  TO  680 
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C0-118 
Cl— 478 
C2-912 
GO  TO  680 
C0- 2676675 
Cl— 5636789 
C2- 2968036 
GO  TO  680 
C0— 141471 
Cl -308123 
C2— 166107 
GO  TO  680 

DHPC- < C0+ ( C1*AM I C ) + ( C2* ( AM I C**2 ) ) ) *THETA 

HC( I )-ALTC I )+DHPC 

DPPPS- ( 3 . 61382E-5*DHPC ) /THETA 

DMPC- ( ( 1+ ( . 2* ( AM I C**2 ) ) ) *DPPPS ) /C 1 . 4*AM I C ) 

AMCC  I ) -AMIC+DMPC 

UTAS ( 1) - 1 . 6878*AMC ( 1) *38 . 96763XSQRT ( THETA*288 .  15  ) 

CONTINUE 

THIS  PORTION  COMPUTES  TRUE  AOA  BY  THE  ANGLE  METHOD.  FLIGHT 
PATH  ANGLE  IS  CALCULATED  BY! 

FLIGHT  PATH  ANGLE  CFPANG)  =  I NO  SIN  CUERTICAL  UEL  (HDOT)  / 
TRUE  AIRSPEED  CUT AS)] 

TRUE  AOA  IS  CALCULATED  BY'  AOA  -  PITCH  ANGLE  -  FLT  PATH 
AOA  IS  THEN  AUERAGED  OUER  A  TIME  INTERNAL  FOR  USE  AS 
THE  INITIAL  AOA  FOR  THE  ITERATIUE  METHOD 

DO  750  1=1, NA 
I Ft  I  .EQ.  1 )GO  TO  720 
IFCI  . EQ.  NA)GO  TO  720 
DTIME-TI C 1+1 )-TI C 1-1 ) 

DALT-HCCI+D-HCC  J-l) 

HDOT ( I) - DALT/DT I ME 

FPANG ( I ) -ASINCHDOTC I )  /UTASC I ) ) 

FPANG ( I ) -FPANG ( I )*57. 29578 

AOA1 ( I) - C TH ( I) -FPANG ( I) ) /COS ( BA C I ) /57 . 29578 ) 

ERRORC I)-AOAl( I)-AOA( I ) 

WRITEC4, 740) I , TI ( I ) , HCCI) . AMCC  I ) . UTASC I ) , HDOTC I ) . FPANGC I) , 
1AOA1C I ) ,AOAC I ) .ERRORC I ) 

FORMATC IX, I 10.9F10. 3) 

CONTINUE 

THIS  PORTION  AUERAGES  THE  ALTITUDE,  MACH,  PITCH, 

PITCH  RATE,  NZ,  AOA-MEAS,  AOA-CALC,  AND  ERROR 

DO  752  I-2.NA-1 
DTIME-TI  CI+D-TICI-l) 

TALT-TALT+  C  HCC I )*DTIME) 

TMACH-TMACH+CAMCC I )*DTIME) 

TTH-TTH+CTHC I )*DTIME) 

TQ-TQ+CQC I )*DTIME) 

TNZ-TNZ+CZNC I ) *DT I ME ) 

TAOA-TAOA+CAOAC I )*DTIME) 

TAOAl-TAQAl+CAOAlC I )*DTIME> 

TERROR-TERROR+CERRORC I ) *DT I ME ) 


V  % 


ttime-ttime+dtime 

CONTINUE 

AALT - TALT/TT I ME 

AMACH-TMACH/TTIME 

ATH-TTH/TTIME 

AQ-TQ/TTIME 

ANZ-TNZ/TTIME 

AAOA-TAOA/TTIME 

RAOA1 -TAOA1/TTIME 

AERROR - TERROR/TT I ME 

DO  7S5  I-2.NA-1 

XALT-AALT-HC( I ) 

XMACH-AMAOHAMCC I ) 

XTH-ATH-THC I ) 

XQ*=AQ-Q( I ) 

XNZ-ANZ-ZNC I ) 

XAOA-AAOA-AOAC  I ) 

XAOA1-AAOA1-AOA1  ( I ) 

XERROR-AERROR-ERROR ( I ) 

XALT*ABS(XALT) 

XMACH-ABS(XMACH) 

XTH-ABSCXTH) 

XQ-ABS(XQ) 

XNZ-PBS(XNZ) 

XAOA-ABS(XAOA) 

XAOAi-ABS(XAOAl) 

XERROR-ABS ( XERROR ) 

IFCXALT  .GT.  YALT5VALT-XALT 
IF(XMACH  .GT.  YMACH ) YMACH-XMACH 
IFCXTH  .GT.  YTH) YTH-XTH 
IF(XQ  .GT.  YGDYQ-XQ 
IFCXNZ  .GT.  YNZ)YNZ«XNZ 
IF(XAOA  .GT.  YAOA)YAOA=XAOA 
IF(XAOAl  .GT.  YAOADYAOAl-XAOAl 
I F ( XERROR  .GT.  YERROR ) YERROR -XERROR' 

CONTINUE 

WRITE  ( 4 , 758 ) AALT . YALT . AMACH . YMACH . ATH . YTH . AQ , YQ , ANZ . YNZ , 
1 AAOA1 , YAOA1 , AAOA . YAOA . AERROR , YERROR 
FORMAT ( / i *  AUG  ALT  -  ’.F10.3.’  +/-  ’.FB.3,’  AUG  MACH 

1F10.3,  ’  +/-  ' i F6. 3, /, ’  AUG  PITCH  «  ’.F10.3, ’  +/-  ’.F6.3 

2*  AUG  PITCH  RATE  -  ’ ,F10.3, ’  +/-  ’.FG.3,’  AUG  NZ  - 

3F10. 3. ’  +/-  ’ . F6. 3. /» ’  AUG  AO A- 1  -  ’.F10.3.’  +/-  ’.F6.3 

4’  AUG  AOA-M  -  ’ ,F10.3, ’  +/-’,F6.3, ’  AUG  ERROR  -  ’ , 

5F10.3,’  +/-  ’ ,F6.3,//) 

THIS  PORTION  COMPUTES  THE  AOA  OF  THE  REMAINING 
POINTS  BY  COMPUTING  FLIGHT  PATH  ANGLE,  THEN 
USING’  AOA  -  PITCH  ANGLE  -  FLIGHT  PATH  ANGLE 

IFCN  .EQ.  NA)GG  TO  772 
DO  770  I-NA+l.N 
IFCI  .EQ.  N)GO  TO  760 
DTIME»TI  ( I+D-TI C 1-1) 

DALT*HC( 1+1 )-HC( 1-1 ) 

HDOTC I )*DALT/DTIME 

FPANG(  I  )-ASIN(HDOT(  I  )ZUTAS(  I ) ) 


J  ■  <•  .1 I  *"*V*W*W^W1H*^*' 


FPANG< I ) -FPANGC I )*57. 29578 

A0AKI)«(TH(I)-FPANG(I))/C0S<BA(I)/57.29578> 

ERROR! I ) “AOA1 ( 1 )— OOflC I ) 

760  WRITE! 4, 765)  I«TI(I)>HC(I)  <  AMC!  I )  ,  UTASC I )  »HDOT(  I )  >  FPANG!  I )  , 

1AOA1 ( I )  #  AOA ( I ) ,  ERROR ( I ) 

765  FORMAT! IX. 110, 9F10. 3) 

770  CONTINUE 

C 

C  THIS  PORTION  COMPUTES  ANGLE-OF-ATTACK  BY 
C  COMPUTING  ALPHA-DOT,  THEN  ADDING  IT  TO 
C  THE  PREUIOUS  AOA  TO  COMPUTE  A  NEW  AOA 
C 

772  DO  800  1*1,  N 
A0A2( 1 ) =AOA( 1 ) 

A0A3C 1 ) “AAOA1 
IFC I  .EQ.  1 )GO  TO  77S 
AUGZN*  ( ZN ( I  )+ZN(  1-1 )  )/2 
AUGBA«!BA! I )+BA( 1-1 ) )/2 
AUGTH- (TH( I )+TH! 1-1 ) )/2 
AUGQ=(Q( I )+Q( 1-1 ) )^2 
AUGTASKUTAS! I )+UTAS! 1-1) )/2 

AZ= AUGZN-  ( COS  ( AUGBA/57 . 29578 )  *COS  ( AUGTH/S7 . 29578 )  ) 

ADOT( I ) =AUGQ--!  AZ*32. 2*57. 29S78/AUGTAS) 

DTIME-TI! I)~TI! I-l) 

A0A2C I )*AOA( 1-1 )+(ADOT( I )*DTIME) 

A0A3! I ) *A0A3!  I-l)  +  (ADOT(  I  )*DTIME) 

C 

C  THIS  PORT I  OH  CORRECTS  AOA  CALCULATED  AT  CG  FOR 
C  PITCH  RATE’-  AOA  CG  *  AOA  WANE  +  !Q  *  X  /  UTAS) 

C  WERE  X  IS  THE  DISTANCE  BETWEEN  THE  WANE  AND  THE 

C  ACCELEROMETER  LOCATION  ON  THE  AIRCRAFT 

C 

X-17 

I F ( ATYPE C 1 )  .EQ.  ’T’)X=25 
A0AT-A0A3! I  )-CQ(  I  )*X/UTAS(  I  ) ) 

ALPHAK=AOA(  I  VAOAT 
ERRORK I )-A0A2( I )-AOA( I ) 

ERR0R2 ! I ) - AOAT-AOA ( I ) 

ADQTM !  I )  *  ( AOA ( I ) - AOA ( I - 1 ) ) /DT I ME 

775  WRITE! 6, 780) I ,TI ( I ) , A  DOT! I ) , ADOTM! I ) , A0A2! I > , AOA! I ) , ERRORK I ) , 

1A0A3! I ) , AOAT , ERR0R2 ( I ) , ALPHAK 
780  FORMATC1X, 110, 10F10.3) 

800  CONTINUE 

C 

STOP 

ENU 
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PROGRAM  KALOPT 

C 

C  PROGRAM  KALOPT  IS  A  FORTRAN  PROGRAM  DESIGNED  TO  OPTIMIZE 
C  ANGL E-OF-ATT ACK  USING  KALMAN  FILTER  EQUATIONS.  THE 

C  PROGRAM  USES  DATA  PARAMETERS  RECORDED  INFLIGHT  ON  AN 
C  AYD I N-UECTOR  DATA  ACQUISITION  SYSTEM.  KALOPT  USES  TIME  CTI ) . 

C  IAS  CAS),  INDICATED  ALT  (ALT),  MEASURED  AOA  (AOA) >  PITCH 
C:  ANGLE  CTH) ,  PITCH  RATE  (Q),  NORMAL  ACCELERATION  CZN),  SIDESLIP 

C  ANGLE  (AOSS) ,  AND  BANK  ANGLE  (BA)  OBTAINED  FROM  DYNAMICS  EUP 

C  AND  EUS  FILES.  AOAOPT  CORRECTS  FOR  PITOT-STATIC  ERRORS  TO 
C  DETERMINE  PRESSURE  ALTITUDE  (HC),  MACH  NUMBER  (AMC),  AND 
C  TRUE  AIRSPEED  (UTAS).  THE  PROGRAM  CALCULATES  AN  "OPTIMAL"  TRUE 
C  AOA  BY  COMBINING  STATE  ESTIMATION  OF  AOA  AND  PITCH  ANGLE  WITH 

C.  ACTUAL  PITCH  ANGLE  MEASUREMENTS.  THE  PROGRAM  WEIGHS  THE  PITCH 
C  ANGLE  MEASUREMENT  AND  ADDS  THE  WEIGHTED  UALUE  TO  THE  ESTIMATED 
C  AOA  AND  PITCH  ANGLE.  THE  UPDATED  AOA  IS  THE  "OPTIMAL"  TRUE  AOA. 
C  THIS  PROGRAM  CAN  BE  USED  WITH  T-38A  OR  RF-4C  DATA. 

l 

BYTE  PNAME ( 12) ,  ENAMEC 12) ,  D0F(7) ,  ATYPEC5)  ,  ATAILO) 

DIMENSION  TI (300) , AS( 300) . ALT(300) , AOAC300) , TH( 300) , Q(300) , 

1ZN ( 300 ), AOSS (300) , BA ( 300 ) , HC ( 300 ) , AMC ( 300 ) >  UTAS ( 300 ) 

DIMENSION  B(2,2),G(2,2),QA(2,2),P(2,2),AK(2>1),XA(2,1),BT(2,2), 
1U(2,1),XP(2,1),GT(2,2) , QGTC2, 2) , GQGT(2, 2) , PP(2. 2) ,H( 1 , 2) . 

2HT ( 2 , 1 ) >  PPHT (2,1), AKZMXH ( 2 , 1 ) , HPP (1,2), AKHPP (2,2) 

BYTE  FILE1( IS) , OF I LEI ( 17) 

BYTE  F I LE2 ( IS ) , OF I LE2 ( 17 ) 

BYTE  FILE3( IS) , 0FILE3C17) 

BYTE  F I LE4 ( 15 ) . OF I LE4 ( 17 ) 

BYTE  FILE6( IS) , 0FILE6C 17) 

BYTE  FILE7C1S) ,0FILE7(17) 

BYTE  FILE8C15) .0FILE8C17) 


DATA 

FILEl/’E’ 

’S’  ,  * 

!’ ,8*0, ’ . 

,  ’E’ 

,  ’U’  , 

,p’  / 

DATA 

FILE2/’E’ 

’S’  ,  ’ 

i’ ,8*0. ’ . 

,  ’E’ 

,  ’U’  , 

’S’/- 

DATA 

FILE3/’E’ 

’S’  ,  ’ 

,8*0,  ’  . 

,  ’D’ 

, ’A’ , ’T’/ 

DATA 

FILE4/’E’ 

’S’  ,  ’ 

■  ,8*0, ’ . 

,  ’D’ 

,  ’A’  ,’TV 

DATA 

FILE6''’  E' 

’S’  ,  ’ 

»’ ,8*0, ’ . 

,  ’D’ 

, ’A’ , ’T’s 

DATA 

F I LE7/ ' E ’ 

•S’ ,  ’ 

i  ’ .8*0, ’ . 

,  ’D’ 

,  ’A’ , ’ T’ / 

DATA 

F I LE8/ ’ E ’ 

'S’  ,  ’ 

i  ’ ,8*0. ’ . 

,  ’  I  ’ 

,  ’N’ , ’P’/ 

C 

WRITE  (5,100) 

100  FORMAT (’  ENTER  NAME  OF  DYNAMICS  FILE>  NNNNN . EUP ’ ,/) 

ACCEPT  110, (FILE1C I ) , 1=4, 11 ) 

110  FORMAT  (8A1 ) 

OPEN  ( UN I T  = 1 ,  NAME  - F I LE 1 .  TYPE  *  ’  UNKNOWN  ’  ) 

WRITE  (5,120) 

120  FORMAT ( ’  ENTER  NAME  OF  DYNAMICS  FILE!  NNNNN. EUS’ , /) 

ACCEPT  110, (FILE2(I), 1-4,11) 

OPEN ( UN  I T - 2 , NAME = F I LE2 . TYPE  = ’ UNKNOWN ’ ) 

WRITE  (5. 125) 

125  FORMAT (’  ENTER  NAME  OF  INPUT  FILE'  NNNNN. I NP’/) 

ACCEPT  110, (FILE8(1), 1-4, 13 ) 

OPEN (UNIT-8,  NAME =F I LEB , TYPE  = ’ UNKNOWN ’ ) 

WR I TE ( 5 , 13Mj 

130  FORMAT (’  ENTER  NAME  FOR  DATA  FILE  1=  NNNNNN.DAT’ ,/) 

ACCEPT  110,  (FILE3(I),  1-4,11) 

0PEN(UNIT=3, NAME-FILE3, TYPE- ’  USKNQWN’ ) 
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WRITE(5» 140) 

140  FORMATC  ENTER  NAME  FOR  DATA  FILE  2'  NNNNNN.DAT’ »/) 

ACCEPT  110, (FILE4( I), 1-4,11) 

OPEN(UN I  T=4 « NAME-F I LE4 , TYPE” ’ UNKNOWN ’ ) 

WRITE(5, 150) 

150  FORMAT ( ’  ENTER  NAME  FOR  DATA  FILE  3«  NNNNNN.DAT’ ./) 

ACCEPT  110, (FILE6(1), 1-4,11) 

OPEN ( UN I T - 6 , NAME - F I LE6 , TYPE = ’ UNKNOWN ’ ) 

WRITE(5, 155) 

155  FORMAT  ( *  ENTER  NAME  FOR  DATA  FILE  4-*  NNNNNN.DAT’,/) 

ACCEPT  110, (FILE7C1), 1-4.11) 

0PEN(UNIT-7,NAME*FILE7, TYPE- ’UNKNOWN’ ) 

WRITE(5. 160) 

160  FORMAT  (’  ENTER  NUMBER  OF  DATA  POINTS  DESIRED  ’  ,/) 

READ(5, 170)  N 
170  FORMAT  (IS) 

WRITE(5. 100) 

180  FORMAT  (’  ENTER  FIRST  LINE  OF  DATA  DESIRED  ’  ,/) 

READ (5, 170)  M 
WRITE(5, 190) 

190  FORMAT ( ’  ENTER  NUMBER  OF  INITIAL  CONSTANT  AOA  LINES  ’  ,/) 

READC5, 170)NA 
WRITE(5,200) 

200  FORMAT ( ’  ENTER  CORRECTIONS  FOR  PITCH  RATE  AND  NZ<  XXX.YYY’ ,/. 

1’  CORR  ARE  FROM  GROUND  BLOCK,  ABOUE/BELOW  0  FOR  Q’ ,/, 

2’  ABOUE/BELOW  1  FOR  NZ’,/) 

READ ( 5 , 2 1 0 ) QCORR , ZNCORR 
210  FORMAT (2F10. 3) 

C 

C  THIS  PORTION  READS  DATA  FROM  THE  DYNAMICS  EUP  AND  EUS  FILES 
C  AND  FORMATS  FOUR  DATA  FILES  (XXXX.DAT)  WHERE  THE  PROGRAM 
C  RESULTS  WILL  BE  SENT.  DATA  FILE  #1  RECORDS  THE  RAW  DATA  FROM 
C  THE  EUP  AND  EUS  FILES.  DATA  FILE  #2  RECORDS  THE  OUTPUT  OF 
C  THE  KALMAN  FILTER  PROGRAM <  "OPTIMAL"  UALUES  FOR  AOA  AND  PITCH 
C  ANGLE,  AND  THE  ERROR  BETWEEN  THE  "OPTIMAL"  AOA  AND  MEASURED 
C  AOA.  DATA  FILE  #3  RECORDS  THE  UALUES  OF  THE  PCT)-  AND  KCT) 

C  MATRICES.  DATA  FILE  #4  RECORDS  THE  UALUES  OF  THE  P(T)+ 

C  MATRIX  AND  THE  TRUE  RMS  ERROR. 

C 

READ  ( 1 , 220 ) PNAME , ENAME , DOF 
220  FORMAT  (8X, 12A1, 15X, 12A1 , 21X, 7A1) 

WRITE  ( 3 , 230 ) PNAME , ENAME , DOF 
WRITE  (4. 230)PNAME, ENAME, DOF 
WRITE  (6, 230) PNAME,  ENAME,  DOF 
WRITE  (7.230)PNAME, ENAME, DOF 

230  FORMAT  (’  PILOT < ’ > IX, 12A1.5X, ’ ENGINEER < ’ » IX, 12A1 ,5X, 
l’DATE  OF  FLIGHT « ’ , IX, 7A1 ) 

READ  ( 1 , 240 ) ATYPE , ATA I L 
240  FORMAT  ( 11X, 5A1 . 13X, 3A1 ) 

WRITE  (3, 250) ATYPE, ATA I L 
WRITE  ( 4 , 2S0 ) ATYPE , ATA I L 
WRITE  (6,2S0)ATYPE> ATAIL 
WRITE  (7, 250) ATYPE, ATAIL 

250  FORMAT ( ’  A/C  TYPE' ’, IX, 5A1.5X, ' TAIL  #>’,1X,3A1) 

WRITE  (3,260) 

260  FORMAT  (//,7X» ’LINE’ »6X, ’TIME’ «2X, ’AIRSPEED' ,2X, 'ALTITUDE*  *7X, 
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nnnnnn, 


i'AOA'  >5X* 'PITCH' «2X, 'PITCH  RT* ,8X, ’NZ’ ,6X, ’AOSS’ ,2X, ’ROLL  PNG’) 
WRITE  (4,270) 

270  FORMAT  C//.7X,  'LINE'  ,6X,  ’TIME’, IX,  ’PRESS  ALT’  ,6X,  ’MACH*  ,3X,  ’TRUE 
i  AS’  ,5X,  ’PITCH’  ,3X,  'PITCH  +’,5X,’A0A  +  ’  ,5X,  ’Q*X/V’  ,2X,  ’  AOA-TRUE’ 
22X,*A0A-UANE’,5X, ’ERROR ’.3X,’K-ALPHA') 

WRITE  (6,280) 

280  FORMAT  (//, 7X, ‘LINE’ , 6X, ’ TIME’ , 5X. ' P-(l) ' ,SX, ' P-(2) ’, 5X, *P-(3) * . 
15X, *P— (4) *  »6X, ’K(l)’ ,6X, ’K(2)’  ) 

WRITE  (7,285) 

285  FORMAT  (//,7X, 'LINE’ , 6X, ’ TIME’ ,5X, ’ P+( 1) ’ ,SX, ’ P+(2) ’ ,5X, ’ P+(3) ’ , 
15X, ’P+(4)’ ,4X, ’ZABOUE’ ,4X, ’ZBELOW’ ,4X, ’ZA-TOT’ ,4X, ’ZB-TOT’ ,6X, 
2’DIFF’ ,7X, ’SUM’ ,5X, ’PSQRT’ ) 

WRITE  (3,290) 

290  FORMAT  (16X. ’ (SEC) ’, 3X, ’ (KNOTS) ’, 4X, ’ (FEET) ’ ,5X, ’ (DEG) ’ , 

15X, ’ (DEG)’ , IX, * (DEG/SEC)’ ,7X, ’ (G) ’ , 5X, ’ (DEG)’ ,5X, ’ (DEG)’ ,/) 

WRITE  (4,300) 

300  FORMAT  (16X, ’(SEC)', 4X, ’(FEET)’, 12X, ’(FT/SEC)’, 5X, 

1’ (DEG) ’ ,5X, ' (DEG) ' , 5X, ’ (DEG) ’ ,5X, ’ (DEG) ’ ,SX, ’ (DEG) ’ , 

25X, ’ (DEG)’ ,5X, ’ (DEG)’ ,/) 

WRITE(6,310) 

310  FORMAT  (16X, ’ (SEC)’ ,/) 

WRITE(7, 310) 

READ  (1,320) 

320  FORMAT  (/./,✓) 

READ  (2,330) 

330  FORMAT  (/,/.✓,✓,/) 

DO  350  J-l.M 
READ( 1,340) 

READ(2, 340) 

340  FORMAT  (  ) 

350  CONTINUE 

DO  390  I-l.H 

READ(1,360)TI(I),AS(I).ALT(I),AOA(I),TH(I),Q(I).ZN(I) 

360  FORMAT  (31X,F10.3,F10.3,F10.3,F10.3,F10.3,F10.3,20X,F10.3) 
READ(2.370)AOSS( I ) , BA( I ) 

370  FORMAT  (31X.F10.3, 10X.F10.3) 

C 

C  THE  NEXT  TWO  STEPS  CORRECT  PITCH  RATE  AND  NORMAL  G  FOR  BIAS 
C  FOUND  ON  THE  GROUND  (BIAS  COMPUTED  FROM  0  DEG/SEC  FOR  PITCH 
C  RATE  AND  1  G  FOR  NORMAL  G) 

C 

Q( I )-Q( I )~GCORR 
ZN( I ) *ZN( I )-ZNCORR 

WRITE(3, 380) I,TI(I).AS(I ) , ALT( I ) , AOA( I ) ,TH( I ) ,Q( I ) ,ZN( I ) , AOSS( I ) 
1BA( I ) 

380  FORiAT( IX, I 10.9F10. 3) 

390  CONTINUE 

THIS  PORTION  PERFORMS  THE  PITOT-STATIC  CORRECTIONS  TO  COMPUTE 
PRESSURE  ALT,  MACH  NUMBER,  AND  TRUE  AIRSPEED.  IT  USES  THE 
STANDARD  PITOT-STATIC  EQUATIONS  AND  ERROR  COEFFICIENTS  FROM  THE 
USAF  TEST  PILOT  SCHOOL  FILES. 

DO  710  I-l.N 

IF(ALTd)  .GT.  36089) GO  TO  500 
DELTA- ( l-(6. B7S59E-6*ALT( I )) )X»5. 2561 


THETA-l-<6. 87S59E-6*ALT( I ) ) 

GO  TO  G10 

500  DELTA - .  22337* EXP < -4 . 80635E-5*  C  ALT ( I ) -36089 ) ) 

THETA- . 751674 

510  IFCAS(I)  .GT.  661. 48) GO  TO  520 

QCICPA- ( ( l+( .2*( (AS( I)/661 . 48)**2) ) )**3. 5)-l 
GO  TO  530 

520  QCICPA«((166.922*(<AS<I)/'661.48)**7))/(((?*((AS(I)/661.48)**2)> 

1-1 )**?  A) )-1 

530  QC I CPS-QC I CPA/DELTA 

AMIC-SQRT(S*(  ( (QCICPS+1)**(  .2857143)  )-i) ) 

IF(ATYPE(1)  .EQ.  *T' ) GO  TO  600 

C 

C  POSITION  ERROR  COEFFICIENTS  COME  FROM  THE 
C  USAF  TEST  PILOT  SCHOOL  PITOT-STATIC  FILES 
C 

DHPC-(~90?+d0270*AMIC)+(-4449S*(AMIC**2))+(93931*(AMIC**3)) 

1 + ( -95540*  C  AM I C**4 ) )  +  ( 37208* ( AM I C**5 ) ) )*THETA 
I F  ( AM  I C  .GT.  DDHPC-0 
GO  TO  700 

600  IF(AMIC  .LT.  . 95S)GO  TO  620 
IFCAMIC  .LT.  . 967)G0  TO  640 
IFCAMIC  .LT.  1.02S)GO  TO  660 
C0=- 46325 
Cl -82503 
C2— 36667 
GO  TO  680 
620  C0-118 

Cl —478 
C2-912 
GO  TO  680 
640  C0- 2676675 

Cl— 5636789 
C2* 2968036 
GO  TO  680 
660  C0— 141471 

Cl -308123 
C2— 166107 
GO  TO  680 

680  DHPC- CC0+CC1*AMIC)+CC2*CAMIC**2) ) )*THETA 
700  HC ( I ) -ALT ( I ) +DHPC 

DPPPS- (3. 6 1 382E-5*DHPC )  /THETA 

DMPC-C  Cl+C .2*CAMIC**2) ) )*DPPPS)/C1 . 4*AMIC) 

AMC C I ) = AM I C+DMPC 

UTAS C I ) » 1 . 6878*AMC (I ) *38 . 96763*SQRT ( THETA*288 . IS ) 

710  CONTINUE 

C 

C  THIS  PORTION  COMPUTES  TRUE  AOA  BY  THE  RELATIONSHIP)  TRUE  AOA  - 
C  PITCH  ANGLE  -  FLIGHT  PATH  ANGLE.  FLIGHT  PATH  ANGLE  IS 
C  COMPUTED  BY)  FLIGHT  PATH  ANGLE  (FPANG)  -  I NO  SIN  tUERTICAL 
C  UEL  (HDOT)  /  TRUE  AIRSPEED  COTAS) ].  THE  TRUE  AOA  IS  THEN 
C  AUERAGED  OUER  A  TIME  INTERNAL  TO  USE  AS  AOAC0)  IN  THE  STATE 
C  ESTIMATION.  THE  DATA  USED  HERE  IS  UNFILTERED. 

C 

IFCNA  .EQ.  0)AAOA-AOA( 1 ) 

IFCNA  .EQ.  0)GO  TO  790 
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DO  750  I-2.NA-1 
DTIME*T!  ( I+D-TI  ( 1-1) 

DALT-HCCl+D-HCCI-l) 

HDOT-DALT/DTIME 
FPANG-ASINCHDOT/VTASC  1 ) ) 

FPANG - FPANGXS7 .29578 
AOAl-THC I ) -FPANG 
C 

C  THIS  PORTION  WERAGES  THE  CALCULATED  AOA  TO 
C  USE  AS  AOAC0)  IN  THE  STATE  ESTIMATION 
C 

TTIME-TTIME+DTIME 
TAOA-TAOA+ (AOAl*DT I ME ) 

750  CONTINUE 

AAOA-TAOA/TTIME 

C 

C  THIS  SECTION  PRINTS  OUT  THE  INITIAL  LINE  OF  DATA 

C 

790  URITE(4,800)TIC1).HC(1),AMC(1),UTASC1),TH<1),TH(1).AAOA.AOA(1) 

800  FORMATC10X, ’1’ ,7F10.3,20X,F10.3) 

UIRITE(6,810)Tl  ( 1 ) 

810  FORMATC10X. ’1’ .F10.3) 

URITE(7,B10)TI(1) 

C 

C  THIS  PORTION  READS  IN  THE  UALUES  FOR  THE  MATRICES  PCT0).  QCT) , 

C  AND  RCT)  [CALLED  P.QA.AND  RI  FROM  AN  INPUT  FILE  CALLED  XXXX. INP 

C 

DO  910  1*1,2 

READ (8 >900) PC I,1)«PCI,2) 

900  FORMATC11X,F10.3,10X,F10.3) 

910  CONTINUE 

DO  930  1*1,2 

READ(8,920)QA( 1 , 1) ,QA( 1 ,2) 

920  FORMAT(11X,F10.3« 10X.F10.3) 

930  CONTINUE 

READC8,940)R 
940  FORMAT (11X,F10.3) 

C 

C  THIS  PORTION  INITIALIZES  THE  STATE  MATRIX  XCT0)  AND  THE 
C  MATRICES  BCT),  GCT),  H(T) ,  ANDUCT). 

C 

XACl.D-AAOA 

XAC2.1)»TH<1) 

DO  2000  I-2.N 
T=TICI)-TICI-1) 

BC1, 15-1 

B ( 1 . 2 ) = ( -32 . 2 ) *57 . 29S78/UTAS (I ) 

BC2, 1 )  =1 
G(  1 , 1 )  =  1 

SCI, 2)*=(-32. 2)^57. 29578/UTAS(  I ) 

GC2, 1 )  =1 
HC1,2) =1 
U(l> 1)“Q(  I ) 

U ( 2 . 1 ) - ZN ( I) -COS ( TH ( I ) /57 . 29578 ) 

C 

C  THIS  PORTION  PROPOGATES  THE  STATE  MATRIX  FROM  THE  LAST  UPDATE 


C  X(T+)  TO  THE  POINT  PRIOR  TO  THIS  UPDATE  X(T-) 

C 

DO  1020  J-1,2 
DO  1010  K-1.2 
BTCJ.K)»B(J,K)*T 
1010  CONTINUE 
1020  CONTINUE 

DO  1030  J-1,2 

XP(J,l)-XA(J.l)+(BTCJ.l)*U(l,l))+(BT(J,2)*U(2,i:>) 

1030  CONTINUE 

C 

C  THIS  PORTION  PROPOGATES  THE  COUARIANCE  MATRIX  PIT)  FROM  THE 
C  LAST  UPDATE  P(T+)  TO  THE  POINT  PRIOR  TO  THIS  UPDATE,  P(T-) 

C 

DO  1050  J-1,2 
DO  1040  K-1,2 
GT(K,J)«G(J,K) 

1040  CONTINUE 
1050  CONTINUE 

DO  1070  J-1,2 
DO  1060  K-1,2 

GK3T(K.J)-CQA(K.l)*GTa.J))+(QACK,2)*GT(2,J)) 

1060  CONTINUE 
1070  CONTINUE 

DO  1090  J-1,2 
DO  1080  K-1,2 

GQGT(J,K)=(G(J,1)*QGT(1,K))+CG<J.2)*QGT(2.K>) 

1080  CONTINUE 
1090  CONTINUE 

DO  1150  J-1,2 
DO  1140  K-1,2 

PP(J,K)»P(J,K)+<T*GQGT(J,K)) 

1140  CONTINUE 
1150  CONTINUE 
C 

C  THIS  PORTION  COMPUTES  THE  GAIN  MATRIX,  KCT) 

C 

DO  1210  J-1,2 
HT(J,1)-H(1, J) 

1210  CONTINUE 

DO  1230  J-1,2 

PPHT(J,l)-(PP<J,l)*HT(l.l))+CPPCJ,2)*HTC2,i)) 

1230  CONTINUE 

HPPHT- CHC 1 . 1 )*PPHTC 1 , 1 ) )+CH< 1 , 2)*PPHT(2, 1 ) ) 

HPPHTR-HPPHT+R 
DO  1270  J-1,2 
AK( J, 1 ) -PPHTCJ, 1 )/HPPHTR 
1270  CONTINUE 
C 

C  THIS  PORTION  PERFORMS  THE  UPDATE  OF  THE  STATE  MATRIX,  X(T). 

C  TO  TAKE  IT  FROM  X(T-)  TO  X(T+)  BY  WEIGHING  THE  PITCH  ANGLE 
C  MEASUREMENT  (Z)  AND  ADDING  IT  TO  X(T-) 

C 

Z-THCI) 

HXP«<H<1,1)*XPC1.1)>+<H(1,2)*XPC2.1)) 

ZMXHP-Z-HXP 
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DO  1320  J-1.2 

AKZMXHC  J,  1 )  -  J,  1  )*ZMXHP) 

XA(J.l)-XP(J.l)+AKZMXH(J.l) 

1320  CONTINUE 
C 

C  THIS  PORTION  UPDATES  PCT)  FROM  PCT-)  TO  PCT+) 

C 

DO  1400  K-1,2 

HPP(1#K)-(H(1. 1)*PP(1«K) )+(H(1.2)#PP(2«K) ) 

1400  CONTINUE 

DO  1430  J-1.2 
DO  1420  K-1.2 

AKHPP(J,K)-(AK(J,1)*HPP(1,K)) 

P( J. K) -PP( J. K)-AKHPP( J. K) 

1420  CONTINUE 
1430  CONTINUE 

C 

C  THIS  PORTION  CORRECTS  AOO  CALCULATED  AT  CG  FOR 
C  PITCH  RATE'  AOA  CG  -  OOA  UANE  +  (Q  *  X  /  UTOS) 

C  WHERE  X  IS  THE  DISTANCE  BETWEEN  THE  UANE  AND  THE 

C  ACCELEROMETER  LOCATION  ON  THE  AIRCRAFT 

C 

X-17 

IFCATYPE(l)  . EQ.  ’T’)X-2S 
AOAT-XAC 1 . 1  )-CQ(  I  )*X/UTAS(  I  )  ) 

QXU-Q(I)*X/UTAS(I) 

C 

C  THIS  PORTION  COMPUTES  THE  SQUARE  OF  THE  ERROR  BETWEEN 
C  PITCH  ANGLE-CALC  AND  PITCH  ANGLE-MEAS  TO  COMPUTE  THE 
C  TRUE  RMS  ERROR 
C 

ZA-XAC2. 1)-TH( I ) 

IF(ZA  .LT.  0)GO  TO  1500 

ZA-Z0**2 

ZB-0 

ZAT-ZAT+ZA 
GO  TO  1510 
1500  ZB-ZA**2 
ZO-0 

ZBT-ZBT+ZB 
1510  DIFF-ZAT-ZBT 

SUM-(ZA+ZB)**.5 

PSQRT-P(2,2)**.5 

C 

C  THIS  PORTION  COMPUTES  THE  ERROR  BETWEEN  THE  “OPTIMAL"  TRUE  AOA 
C  AND  THE  MEASURED  AOA  AND  THEN  PRINTS  OUT  THE  RESULTS  OF 
C  THE  KALMAN  FILTER  PROGRAM. 

C 

ERROR- AOAT-AOAC I ) 

ALPHAK-AOAC I )/AOAT 

WRITE! 4. 1800) I .Tl ( I ) . HCC I ) . AMCC  I ) .UTASC I ) .THC I ) > 

1XA ( 2 . 1 ) , XA (1 . 1) . QXU , AOAT . AOA ( I) , ERROR . ALPHAK 
1800  F0RMATC1X. I10.12F10.3) 

WRITEC6. 1010) I.TI(I).PP(l.l)iPP(li2).PP(2.1).PP(2f2)» 

1AKC 1.1). AKC2. 1 ) 

1010  F0RMATC1X. I10.7F10.3) 


URITE(7,1820)I,TI(I),P(l.l).Pa,2>,P<2,l),P(2,2) 
1 ZA .  ZB . ZAT . ZBT . D I FF . SUM i PSQRT 
1820  FORMATdX.  U0.12F10.3) 

2000  CONTINUE 
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Block  11.  (Cont)  ANGLE-OF-ATTACK  POSITION  ERROR  FROM  FLIGHT  TEST  DATA 

Block  19.  (Cont)-^- error  testing.  The  method  was  accurate  enough  to  identify 
a  hysteresis  error  in  the  T-38A's  AOA  sensor  of  +/-  0.5  degrees,  which  was 
confirmed  by  ground  calibration. 
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